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Abstract 



Aims. The detection of very hot plasma in the quiescent corona is important for diagnosing heating mechanisms. The presence and 
the amount of such hot plasma is currently debated. The SphinX instrument on-board CORONAS-PHOTON mission is sensitive to 
X-ray emission well above 1 keV and provides the opportunity to detect the hot plasma component. 

Methods. We analyzed the X-ray spectra of the solar corona collected by the SphinX spectrometer in May 2009 (when two active 
regions were present). We modelled the spectrum extracted from the whole Sun over a time window of 17 days in the 1.34 - 7 keV 
energy band by adopting the latest release of the APED database. 

Results. The SphinX broadband spectrum cannot be modelled by a single isothermal component of optically thin plasma and two 
components are necessary. In particular, the high statistics and the accurate calibration of the spectrometer allowed us to detect a 
very hot component at ~ 7 million K with an emission measure of ~ 2.7 x 10 44 cm" 3 . The X-ray emission from the hot plasma 
dominates the solar X-ray spectrum above 4 keV. We checked that this hot component is invariably present both at high and low 
emission regimes, i.e. even excluding resolvable microfiares. We also present and discuss a possible non-thermal origin (compatible 
with a weak contribution from thick-target bremsstrahlung) for this hard emission component. 

Conclusions. Our results support the nanoflare scenario and might confirm that a minor flaring activity is ever-present in the quiescent 
corona, as also inferred for the coronae of other stars. 

Key words. Sun: corona - Methods: observational - Techniques: spectroscopic 
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1. Introduction 

The solar corona is X-ray bright. The bulk of the emission is 
thermal, with a spectrum made of a multitude of emission lines, 
e.g. of iron, overlapping a bremsstrahlung continuum, character- 
istic of an optically thin plasma at 1-3 MK. Most of the emission 
comes from plasma confined in closed magnetic tubes anchored 
in the photosphere, the coronal loops. Hotter plasma, to tem- 
peratures well above 10 MK, is typically detected during very 
impulsive flare events, occurring in localized coronal regions. 

Recently, the high sensitivity of new imaging and spec- 
tral instrumentation has allowed to refine the investigation of 
the thermal structure of the corona. One question is the pres- 
ence of minor very hot component at temperature above 6-8 
MK even outside of macroscopic flare events. This question 
is linked to the more basic question of the presence of flar- 
ing activity down to very small scales, to constitute a major 
steady component of the heating extending to the entire so- 
lar corona. Although very faint, such component has been de- 
tected in the recent past in active regions both with filter ra- 
tios or differential emission measure (DEM) reconstru ction from 
imagin g broad-band instru ments (Hinode/X RT, e. g jReale et alj 
2009aflbl: lMcTiernanll2009t ISchmelz et al.ll2009h and with spec- 
troscopic line a n alysis (single line or DEM recon s truction, e. g.. 



Ko et all l2009t iPatsourakos & Klimchukl 12009: Sh estov et alj 



20iri Svlweste r et alJl2010h . More solid confirmation has been 



found very recently (iTesta & Realdl20l2l) . The weight of this 
very hot component has been estimated to be a few percent at 
most of the major c omponent around 3 MK in an active region 
dReale et alj|200 9a). However, there are reasons to believe that 
the amount of very hot plasma may vary from active region to 
active region and it may be smaller in less prominent active re- 
gions (ITesta etalJl20lU e.g.). 

Yet, another question becomes the weight of this hot com- 
ponent over long timescales and over the global coronal budget. 
Some hints come fro m Sun-as-a-Star analysis with Yohkoh data 
( Argi roffi et all 2008) but with large uncertainties in the hotter 
part. More tight constraints are then needed in t he hot tail. 

Th e Solar PHotometer IN X-rays (SphinX) dSvlwester et all 
I2008L iGburek et all 120111) . part of the THESIS package on 
board the CORONAS-PHOTON solar mission, is appropriate 
to investigate this issue: it was a broad-band spectrometer (1- 
15 keV), with moderate spectral resolution but sensitive up 
to hig h energies, enough to detect possible small hot compo- 
nents. [SyiwesteretaLl d2012l) (hereafter SKG12) analyzed coro- 
nal spectra measured with SphinX at the last solar minimum 
and obtained that they are well-described by a thermal spectrum 
at about 1.7-1.9 MK, with significant flux only in the 1-3 keV 
range. 

Here we analyze the full-disk unresolved corona during a 
low solar activity period but in the presence of active regions, 
with the aim to search for the detection of a possible minor hot 
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component. In this analysis, we take advantage of the possibil- 
ity to integrate over several observation days. This allows us to 
achieve a very high signal-to-noise ratio and to access to higher 
energies (> 3 keV), where the steep coronal spectrum yields a 
very low count rate, but which is crucial to be sensitive to a hot 
emission measure component. 

Although we include active regions in the analysis, we ad- 
dress the non-flaring corona and therefore in the data selection 
we exclude any major flare event. We will discuss about the role 
of minor, still resolved, flaring components (microflares). 

Going to such high energies one may wonder about the pos- 
sible role of non-thermal emission. High-energy spectrum com- 
ponen ts are for instance typically detected with RHESSI during 
flares ( Veronig et al. 2006t iKrucker et al. 20081 e.g.). In the case 
of lower activity, we may hope in principle to detect non-thermal 
emission in the SphinX spectral range. So we pay attention also 
to this issue. 

The paper is organized as follows: in Sect. [2] we describe 
the data and the data analysis procedure; in Sect. [3]we show the 
results of the spectral analysis (and of the count-rate resolved 
spectral analysis). Our conclusions are discussed in Sect.|4] 




2. Data Analysis 

The spectra presented in this paper were collected by the Dl de- 
tector of the SphinX spectrometer. The detector energy range 
is 1.2 - 14.9 keV and the spectral resolution is ~ 460 eV. 
F urther details on the instrumental characteristics can be found 
in lGburek etail d201ll) and in SKG12. 

We considered all the observations between 2009 
May 7 and 2009 May 24 (the FITS event files 
are available at the SphinX data catalogue website: 
|http : //156 . 17 . 94 . l/sphinx_ll_c atalogue/SphinX _cat_ 
The choice of this time window was dictated by the relatively 
high solar X-ray flux associated with the presence of active 
regions and by the lack of significant flaring events. The active 
region AR11017 (hereafter AR1) was visible between May 
7 and 22, while a second, dimmer, active region (hereafter 
AR2) was visible from May 15 to the end of our time window. 
Figure Q] shows the Hinode/XRT synoptic images of the Sun 
(Ti-poly filter) in May 7, May 15, and May 23 together with 
the corresponding positions of AR1 and AR2. 

To reduce a possible contamination from solar flares, we 
carefully inspected the light-curves of all the observations and 
excluded observations SPHINX_090509_192730J254950 and 
SPHINX_090514_062135.064952 to remove a B1.0 and a A5.9 
flare, respectively. 

For each event file we extracted the corresponding spec- 
trum by applying a very conservative screening criterium to re- 
move particle-related events, spurious mesurements, and non- 
GTI events. Namely, we selected FLAG=0 events only (see the 
SphinX web site for a detailed description of individual flags). 
We then summed the single spectra to obtain a total spectrum 
with more than 3.9 x 10 7 counts and 57 ks of filtered exposure 
time. 

We adopted the response matrix 

SPHINX_RSP-256_nom_D l .fts and perfo rmed spectral analysis 
by using XSPEC V12.7 (lArnaudfl996T) . We performed gain 
fits within Xspec to account for possible weak variations of 
the response gain as an effect of fluctuations in the detector 
temperature during the 17 days spanned by our observations. 
This procedure shifts the energies E on which the response 
matrix is defined according to the formula E' = (E/g\) - g2, 
where (g\,g2) are free parameters in the fit. For all the models 



.main, 
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Figure 1. Hinode/XRT synoptic images of the solar corona 
(Ti-poly filter) in 2009 May 7 (upper panel), 2009 May 15 (cen- 
tral panel), and 2009 May 23 (lower panel). The position of ac- 
tive regions AR1 and AR2 is marked. 



described below we found g\ ~ 1.025, g2 0.03, in agreement 

with the expectations of small variations in the response matrix. 
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Spectral analysis was performed in the 1.34 - 7 keV en- 
ergy band (i. e. instrument channels 24 - 121) to minimize 
calibration issues (1.2 — 1.34 keV) and contaminations from 
the electronic noise (7 - 14.9 keV). The total spectrum was 
modelled by adopting the APEC spectral code (optically thin 
coron al plasma in collisional ionization equilibrium. ISmith et al.l 
2001) based on the 2.0 release of the AtomDB database (see 
http://www.atomdb.org). Because of the relatively limited 
SphinX spectral resolution it was not possible to clearly iden- 
tify line emission features in the spectra and to perform accurate 
diagnostics of the abundances. Therefore, we s et the abundance 
table in the spectral model to that reported by Feldman ( 1992) 
for the solar upper atmosphere. We added a systematic 5% error 
term in our spectral fittings to account for the estimated uncer- 
tainties in the calibration of the spectrometer. 

3. Results 

We first verified that the total SphinX broad-band spectrum can- 
not be properly described by a single isothermal component. A 
fit of this simple model to the spectrum provides an unacceptable 
reduced^- 2 = 7.06 (with 93 d. o. f.) and significant residuals are 
visible in the high-energy tail of the spectrum, as shown in the 
upper panel of Fig. [2] 

The SphinX spectrum can be well fitted by adding a second 
thermal component (see middle panel of Fig. I2J. The model with 
two components provides a reduced^- 2 = 1.07 (with 91 d. o. f.) 
and a null hypothesis probability > 30%. The 1 - 8 A X-ray flux 
is 1 .4 x 1 0~ 5 erg/s /cm 2 and the 1 - 1 5 ke V luminosity is 4.7 x 1 23 
erg/s. Table Q] shows the best-fit parameters obtained with one 
and two thermal components. We did not obtain significant im- 
provements in the quality of the fits by adding a further thermal 
component. In particular, according to the F-test, the probabil- 
ity that the improvement in the fit is not significant is more than 
three times higher than our threshold at 0.1% and the error bars 
in the best-fit parameters grow larger. We therefore conclude that 
the broadband modeling of the SphinX spectrum requires two 
thermal components, a warm component at ~ 2.7 x 10 6 K and a 
hot one ~ 7 x 10 6 K (hereafter 2T model). 



Table 1. Best-fit parameters obtained by modelling the SphinX 
spectrum with one (IT) and two (2T) isothermal components of 
optically thin plasma, and with one thermal component plus a 
power-law (lT-Pow). All errors are at the 90% confidence level. 
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Figure 2. Upper panel: SphinX spectrum of the solar corona col- 
lected between 2009 May 7 and 2009 May 24 together with its 
best-fit model consisting of one isothermal component and resid- 
uals. Middle panel: same as upper panel with a two-component 
thermal model. The contribution of each component is shown. 
Lower panel: same as middle panel with the lT-Pow model. 



The high-energy tail of the SphinX spectrum can also be fit- 
ted by a non-thermal component. A model with a (warm) ther- 
mal component plus a power-law (lT-Pow model, see lower 
panel of Fig. [2]i provides a x 2 = 94.0 (with 91 d. o. f.) that is 
even lower than that obtained with the 2T model. In this case, 



the values of T\ and EM\ are consistent with those obtained with 
the 2T model, while for the non-thermal component we obtain 
a quite steep photon index F, as shown in Table Q] Therefore, 
with the available data it is not possible to exclude that the hard 
X-ray emitting component is associated with non-thermal pro- 
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cesses and/oris a combination of thermal and non-thermal emis- 
sion. 

SKG12 have analyzed SphinX spectra extracted from the 
lowest-activity periods of the 2009 solar minimum and charac- 
terized by the lack of major active regions. They found that in 
time intervals of exceptionally low X-ray flux, SphinX spectra 
can be modelled by a single thermal component at 1 .7 - 1 .9 x 10 6 
K with emission measure ranging between 4 x 10 47 crrT 3 and 
1.1 x 10 48 cm 3 . We note that in the spectra studied by SKG12 
the 1-15 keV X-ray luminosity is a factor 6-10 lower than that 
in our spectrum. The differences in X-ray luminostity, plasma 
temperature, and spectral model can be related to the presence 
of the active regions AR1 and AR2 in our observations. We 
examined the Hinode/XRT synoptic images of the Sun in the 
thin-Be filter, whose bandwidth is approximatively 0.9 - 2.5 
keV (three observations available in our time window: May 16, 
19, and 20) and found that the contribution (in terms of DN de- 
tected by XRT) of AR1 and AR2 is only -6-9% of the total. 
Therefore, the contribution of the active regions to the total flux 
dominates at high energy only. 

In principle it is possible to ascertain the contribution of AR1 
and AR2 by assuming that the spectrum of the underlying corona 
is described by the model adopted by SKG12. Under this as- 
sumption one can add a third, cooler, component to the 2T model 
(by forcing in the fitting process that its temperature, Tq, and its 
emission measure, EMq, are in the ranges derived by SKG12). 
Unfortunately this approach is not feasible since, as explained 
above, the addition of a new component makes the measures of 
the best-fit parameters less precise. Nevertheless, this 3T model 
yelds T ~ 1.9 X 10 6 K and EM ~ 1 x 10 48 ctrT 3 , while T u T 2 , 
EM\, and EM2 are all consistent with those reported in Table Q] 
at less than one sigma (but with much larger uncertainties). 

Upper panel of Fig. [3] shows the 68%, 90%, and 99% confi- 
dence contour levels of the emssion measure of the warm com- 
ponent, EM\ , versus the emission measure of the hot component 
EM2- The contours indicate that the ratio EM1/EM2 is almost 
constant at ~ 2500. The best-fit value of is entangled with 
EM2, as shown by the 68%, 90%, and 99% confidence contour 
levels in lower panel of Fig. [3] 
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Figure 3. Upper panel: 68%, 90%, and 99% confidence contour 
levels of the emssion measure of the hot component, EM2, ver- 
sus the emission measure of the warm component EM\ obtained 
with the 2T model Lower panel: Same as upper panel for EM2 
and T%. 



3.1. Count-rate resolved spectral analysis 

Figure|4]shows the distribution of the X-ray count-rate observed 
by SphinX in our 17 days time window. The average count-rate 
is of a few hundreds counts s _1 , but a peak at ~ 1400 counts s _1 
is also visible. The peak at high count-rate may be associated 
with microflares not directly visible through a simple inspection 
of the light curve. These microflares may also be responsible for 
the hot component detected in the solar spectrum. To investi- 
gate this issue we performed a count-rate resolved spectral anal- 
ysis by extracting the spectrum of all the events detected at low 
count-rate (< 650 s _1 , hereafter LR spectrum) and the spectrum 
of all the events detected at high count-rate (> 650 s _1 , HR spec- 
trum). Figure [5] shows the LR and HR spectra. 

We found a clear relationship between count-rate and spec- 
tral hardness, with the HR spectrum being significantly harder 
than the LR spectrum. Namely, the soft (1.34-3 keV) to hard 
(3-7 keV) flux ratio is ~ 1000 in the LR spectrum and ~ 500 
in the HR spectrum. This result supports a possible link between 
high count-rate and microflare activity. Nevertheless, though the 
LR spectrum is softer than the HR spectrum, the former cannot 
be modelled by a single thermal component (reduced^ 2 = 3.57 
with 93 d. o. f.) and a proper fit still requires an additional hot 
component (reduced^ 2 = 0.98 with 91 d. o. f.), or a power-law 



'W\.. 



count rate (s" 1 ) 

Figure 4. Distribution of the X-ray count-rate observed by 
SphinX between 2009 May 7 and 2009 May 24. The red vertical 
line indicates the threshold chosen to extract the low count-rate 
spectrum and the high count-rate spectrum (see text). 
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Energy (keV) 

Figure 5. SphinX spectrum of all the events detected at low 
count-rate (< 650 s _1 , black crosses) and at high count-rate 
(> 650 s~' red crosses) between 2009 May 7 and 2009 May 
24. 



(reduced^- 2 = 1.01 with 91 d. o. f.). This clearly proves that the 
hot component cannot be only an effect of microflares, since it 
is present even in regimes of low count-rate, where any possible 
contamination from microflares is removed. 

Table [2] shows the best-fit parameters obtained by fitting the 
LR and HR spectra with one and two thermal components, and 
with one thermal component plus a power-law. Because of the 
lower statistics, error bars are somehow larger than those ob- 
tained for the total spectrum and we still found an entanglement 
between the best-fit values of EM2 and T2 (similar to that shown 
in Fig. [3]). Table [2] shows that the best-fit temperatures obtained 
for the HR spectrum are higher than those obtained for the LR 
spectrum (and from those obtained for the total spectrum, see 
Table[TJ. This is coherent with the HR spectrum showing harder 
X-ray emission. 



4. Discussion and conclusions 

We have illustrated the analysis of full-disk X-ray data collected 
with SphinX spectrometer over a two weeks period when one or 
two active regions were ever present on the solar disk. Although 
SphinX has no spatial resolution, its broad band and sensitivity 
up to high X-ray energies (> 10 keV) make it suited to inves- 
tigate the interesting question of the presence and role of high 
temperature plasma components in the solar corona. 

Although the solar X-ray spectrum is quite steep, integrating 
over two weeks allows us to reach a very high signal-to-noise ra- 
tio well above 3 keV. This makes a real difference from previous 
studies. We analyze such a broad-band spectrum, with one and 
two-thermal component models and we obtain that the spectrum 
is not well-described by a single thermal component and a sec- 
ond hotter component improves the fitting at high significance 
level. In practice, in addition to the large warm component pre- 
viously found in the very quiet sun, that we confirm (T ~ 2.7 
MK), our analysis detects a minor hot component (r»7 MK) 
at extremely high significance. The X-ray emission from the hot 
plasma dominates the solar X-ray spectrum above 4 keV and it 
is clearly associated with the presence of the active regions, that 
were absent in the analysis of SKG12. The temperature of this 
hot component is compatible with that found in other studies of 



Table 2. Best-fit parameters obtained by modelling the LR and 
HR spectrua (see text) with one (IT) and two (2T) isothermal 
components of optically thin plasma, and with one thermal com- 
ponent plus a power-law (lT-Pow). All errors are at the 90% 
confidence level. 
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active regions (Reale et al. 2009a; Testa et al. 201 1). The results 
of present investigations are also in an agreement with the multi- 
temperature analysis (DEM reconstruction) based on the 15 line 
fluxes observed by RESIK spectrometer in the range 3.5 and 6 
A dSvlwester et alJ2010l) : the examination of non flaring periods 
(activity classes from A9 to B5) at the beginning of 2003 indi- 
cated that temperatures of hot DEM component are between 6 
MK and 9 MK and its emission measure is 3 orders of magni- 
tude smaller than that of the cooler component. 

Our fittings show that the emission measure of the hotter 
component is less than 0.1% of the other cooler one. As men- 
tioned in SectQ"! in other active regions there was evidence of a 
hot component of the order of a few percent of the main warm 
one. Although this evidence was subject to important uncertain- 
ties and needed quantitative confirmations, there are reasons that 
may help to understand even this amount of discrepancy. The 
active regions studied previously were quite intense ones. Less 
intense active regions did not show this amount of hot plasma, 
and might be mo re compatible with the amount that we find 
dTesta et al.ll201 ll) . Similarly, here we are close to the solar min- 
imum, and these active regions are not very intense. Also, here 
we are analyzing full-disk observations, including the quiet Sun 
emission. This emission may contribute significantly to the de- 
tected warm emission and much less to the hot one, making the 
ratio hot/warm decrease. Finally, we integrate over two weeks 
and therefore we average spatially but also on a very long time 
baseline, that may include periods when the hot component is 
relatively high and periods when it is not. In other words, we 
provide the time-average of the hot component whereas the other 
measurements were snapshot. This time-averaged value might 
be interesting when evaluating the importance of heating mech- 
anisms that invo lve the systematic presence of hot plasma, as 
nanoflare theory (iParkerll 1 9881; ICargillll 1 994t lKlimchukll2006l) . 

In any case, although at such small amount, the hot plasma 
detected by our analysis is at extremely high significance and 
shows that the data collected with SphinX are highly sensitive to 
this component. 
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One important question about the hot component is whether 
it might be due to minor, still resolvable, microflares or to con- 
tinuous, widespread and very small nanoflares. The Soft X-ray 
Telescope (SXT) on Yohkoh found micr oflares to occur in ac- 
tive regions in the 0.25 - 4 keV band (IShimizulll995l) . Also 
RHESSI observes microflares exclusively in active regions, al- 
though with different rates and in the hard X-rays. Microflares 
observed systematically with RHESSI typically show similar 
spectral features as thos e reported in Fig. [2] a nd Fig [5] although 
in a harder energy band dHannah et al.ll2008l 1201 ll) . The differ- 
ent energy bands correspond to d ifferent fitting param eters, as 
expected from experimental bias (lHannah et al.l l20lH) . We can 
equally attempt a comparison if we either consider microflares 
inc luded in or exclu ded from our data. The microflares analyzed 
by iLin etaT. I (fl98l show nonthermal emission with r = 4 - 6 
above 30 keV, i.e. at considerably higher energies with respect to 
the SphinX energy band. The electron spectrum can be assumed 
to exist above some "cutoff" energy EC , which means the ac- 
celerated electron population appears abruptly at energies where 
the thermal spectrum is negligible. For microflares, the observed 
break tends down to very low energies where there are multiple 
emission lines in the thermal spectra. This makes it very diffi- 
cult to determine EC accurately and to determine the energy in 
a microflare's nonthermal electrons, and this holds also for our 
analysis. A systematic analysis of a large sample of microflares 
observed with RHESSI has shown an average po wer law index 
r - 6 .9, with most events in the range 4 < T < 10 (Hannah et al. 
120081) . Our nonthermal fitting falls in the right tail of this dis- 
tribution, i.e. in the low-energy tail. The microflare spectra col- 
lected by RHESSI are often so steep that they are diffic ult to 
distinguish from the thermal emission (Han nah et al.l2 008). as it 
occurs in our analysis. Thermal fitting yield temperatures abov e 
10 MK with an average of about 13 MK (lHannah et al. 2008). 
These temperatures are significantly higher than those obtained 
for the hot component in our analysis. Part of the dicrepancy 
might be explained with the different sensitivity of the instru- 
ments, due to the different spectral band, that would make them 
give different weights to a broad emiss ion measure distribution 
(Pere s et al. I I27M lHannah et alj I20TTI) . The average emission 
measure obtained from thermal fitting of RHESSI microflares 
is about 2 - 3 x 10 46 cirr 3 , much larger than the emission mea- 
sure of the hot components that we obtain. We should consider 
that our value is obtained as an average value over a long time 
baseline. Since the microflaring is not continuous, for a sound 
comparison we should scale the values to i nclude the microflar e 
duty cycle. If we assume a 10% duty cycle dHannah et al.ll20lH) . 
our value still remains quite lower than expected if microflares 
are included. 

Nevertheless, we have tried to disentangle a possible contri- 
bution from microflare activity with an independent approach, 
namely we analyzed separately the spectra obtained from the 
low and high count-rate periods. We obtained a very robust an- 
swer: the hot component is confirmed at high significance in both 
regimes, with the only difference of noticeably different temper- 
ature values and emission measures, but that do not affect the 
general result. Therefore, we find this hot component at any level 
of emission intensity. 

Finally, we have also examined the possibility that this 
hot component could be instead due to non-thermal emission. 
Although on the basis of other previous studies the expecta- 
tion that the emission is thermal is solid, the data do not al- 
low us to exclude that the excess in the hard energy tail is 
due to non-thermal emission. Fitting with a non-thermal compo- 
nent provides sound results, i.e. a power law with quite a steep 



spectral index, still compati ble with a weak con tribution from 
thick-target bremsstrahlung (IVeronig et al.ll2006l) . Nevertheless, 
the lT-Pow framework would imply the lack of hot X-ray emit- 
ting plasma from active regions at odd with previous studies. 
Discriminating between thermal and non-thermal emission re- 
quires at the same time high sensitivity and good spectral reso- 
lution, to resolve possible high-energy lines such as the Fe XXV 
6.7 keV line. Even in that case we probably could not exclude a 
simultaneous presence of thermal and non-thermal spectral com- 
ponents. We point out that if part of the hard X-ray emission that 
we detected has a non-thermal origin, the EM values for the hot 
component reported in Table Q] and in Table [2] might be consid- 
ered as upper limits. In both cases of thermal and non-thermal 
emission this invariably brings support to the presence of steady 
impulsive heating in active regions. 

In conclusion, this analysis of SphinX data confirms the pres- 
ence of a minor yet steady very hot component in coronal ac- 
tive regions, therefore supporting an important role of impulsive 
processes in bright plasma heating. Future higher resolution and 
high sensitivity X-ray spectroscopy might shed more light on the 
details and nature of this hot emission. 
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